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Abstract: 

A  rotating  ensemble  of  bodies  of  arbitrary  shape  with  angular 
periodicity  scatters  an  electromagnetic  wave  to  produce  a  spectrum 
of  frequency  components  characteristic  of  the  structure  and  its  rate 
of  rotation.  The  spectrum  and  its  properties  are  predicted  through 
electromagnetic  field  theory. 

The  theory  has  been  developed  such  as  to  exploit  the  angular 
periodicity  and  thereby  reduce  the  computational  load  by  a  large 
factor.  A  consequence  of  the  approach  is  that  the  spectrum  is  found 
directly. 

Many  of  the  predictions  have  been  confirmed  by  direct  computation 
of  the  scattered  field  for  a  series  of  rotational  positions  to  simulate 
a  time  series,  followed  by  discrete  Fourier  transforming  to  produce 
the  spectrum. 


1. 


Introduction 


The  intake  of  the  jet  engines  of  an  aircraft  makes  a  major 
contribution  to  its  radar  cross-section.  The  rotation  of  the  engine 
compressor,  each  stage  of  which  exhibits  angular  periodicity,  gives 
this  contribution  a  distinctive  spectrum  by  means  of  which  one 
might  hope  to  identify  the  target  [2,3,4].  In  a  previous  study  [5] 
the  compressor  stage  was  crudely  modelled  as  a  coplanar 
ensemble  of  wires  radiating  from  an  axis.  In  the  present  study 
this  is  generalised  by  replacing  each  wire  by  a  conducting  body  of 
arbitrary  shape.  The  ensemble  of  bodies  exhibits  the  same 
angular  perodicity  as  did  the  wires.  If  each  body  is  given  the 
shape  of  a  compressor  blade  a  stage  of  the  compressor  may  be 
simulated.  The  central  hub  or  axle  from  which  the  blades  radiate 
may  be  included  without  difficulty. 

In  the  previous  study  a  number  of  conclusions  were  drawn  about 
the  characteristics  of  the  spectrum.  These  were  based  on  the  wire 
model,  although  in  most  cases  it  was  intuitively  clear  that  they 
would  be  valid  for  the  more  general  case.  In  the  present  study 
these  intuitions  are  established  rigorously. 

This  study,  like  the  previous,  had  two  principal  goals.  The  first 
was  to  find  a  means  to  calculate  the  field  which  exploited  the 
periodicity  of  the  ensemble  and  thereby  to  reduce  the 
computation  load.  The  second  was  to  express  the  result  as  a 
frequency  spectrum.  In  the  event,  the  means  by  which  the  first 
goal  was  met  automatically  met  the  second. 

The  incident  electric  field  is  expressed  as  a  summation  of 
harmonics  of  the  azimuthal  angle  (measured  around  the  axis). 

This  approach  exploits  the  angular  periodicity  of  the  ensemble  of 
bodies.  Each  harmonic  of  the  field  excites  a  harmonic  component 
of  the  current.  For  each  harmonic  it  is  found  necessary  to 
calculate  the  current  on  only  one  body,  for  only  one  azimuthal 
angle  of  incidence;  the  currents  on  the  other  bodies,  for  other 
directions  of  incidence,  are  found  simply  from  the  first.  In 
consequence,  the  computation  needed  to  solve  for  the  scattered 
field  is  reduced  (compared  to  a  direct  calculation  which  does  not 
exploit  the  angular  periodicity)  approximately  in  proportion  to  the 
number  of  bodies. 

A  further  beneficial  consequence  of  the  harmonic  decomposition 
of  the  incident  field  is  that  the  scattered  field  is  found  directly  as 
a  frequency  spectrum.  The  components  are  lines  at  intervals  of 
the  rotation  rate  multiplied  by  the  number  of  bodies.  This  is 
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consistent  with  the  fact  that  the  ensemble  re-presents  the  same 
aspect  to  the  radar  at  this  rate.  The  spectrum  extends  with 
significant  strength  above  and  below  the  illumination  frequency 
by  an  amount  equal  to  the  Doppler  shift  associated  with  the 
maximum  linear  velocity  of  the  bodies.  Beyond  these  limits  the 
spectral  lines  continue,  becoming  weaker  approximately 
exponentially. 

Many  features  of  the  spectrum  of  the  scattered  field,  including 
those  mentioned  in  the  foregoing  paragraph,  may  be  deduced 
from  the  form  of  the  equations  without  proceeding  to  their 
solution.  At  this  point  in  the  study,  exact  solutions  have  not  been 
calculated.  It  is  to  be  expected  that  further  features  useful  for 
target  identification  may  be  exposed  by  the  exact  numerical 
solutions.  To  this  end  an  approach  to  numerical  solution  has  been 
examined  and  described  in  detail. 

It  is  possible  to  make  predictions  concerning  the  effects  of 
conditions  outside  the  present  simulation,  such  as  the  presence  af 
the  engine  cowling,  the  presence  of  multiple  stages  and  the 
presence  of  stator  blades  in  the  engine.  The  fact  that  the 
compressor  stage  is  surrounded  by  a  cowling  and  not  in  free  space 
implies  that  the  incident  field  is  not  a  plane  wave:  however  the 
true  field  may  be  decomposed  into  a  spectrum  of  plane  waves, 
each  of  which  interacts  with  the  scatterer  in  the  manner  described 
by  the  present  theory.  The  presence  of  multiple  stages  has  the 
result  that  the  scattered  field  has  a  spectrum  which  has  some  of 
the  character  of  a  superposition  of  the  spectra  of  the  separate 
stages,  with  weaker  components  at  intervals  of  the  shaft  rate.  The 
presence  of  stator  blades  does  not  alter  the  frequencies  of  the 
spectrum,  but,  by  introducing  new  mechanisms  of  electromagnetic 
interaction  between  moving  and  stationary  bodies,  causes  the 
spectrum  to  have  significant  strength  at  frequencies  beyond  the 
Doppler  limits  mentioned  above. 

Many  of  the  predictions  and  expectations  arising  from  the  theory 
have  been  confirmed  through  numerical  simulation.  Calculations 
for  simple  shapes  have  been  made  using  the  NEC  program  [6]  for  a 
series  of  angular  positions  and  the  results  submitted  to  a  fast 
Fourier  transform.  Confirmation  has  also  been  obtained  through 
measurement  of  some  wire  structures. 
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2.  The  ensemble  of  bodies 


In  this  section  is  described  the  geometry  of  the  ensemble  of 
bodies,  together  with  the  coordinate  systems  and  some  special 
parameters  to  be  used  below. 

I  conducting  bodies,  with  surfaces  Si,  i=0,l,...,I-l,  are  arranged 
about  an  axis  (in  the  coordinate  system  to  be  used,  the  z  axis)  at 

2,7C 

angular  spacings  of  — ,  see  figure  2.1.  Their  shapes  are  the  same 

and  oriented  such  that  Si  is  generated  by  rotating  So  about  the  z 
axis  throught  the  azimutal  angle 

(D.=i^.  (2.1) 

The  objects  may  be  detached,  as  shown,  or  attached  to  one 
another.  In  the  latter  case  the  individual  bodies  must  be  defined 
by  an  arbitrary  cut  through  the  common  region. 

Rectangular,  cylindrical  and  spherical  coordinates  will  be  used  as 
convenient;  the  system  in  use  will  usually  be  apparent  from  the 
conventional  symbols:  (x,y,z),  (p,(l>,z),  (r,0,(l)). 

If  a  point  on  surface  Sq  has  position  vector  and  coordinates 


with  the  conventional  relations 

Po  =  ‘'‘o  ’  f 0  =  ‘=“‘‘  ?  “ 

Pq  PQ 

I 

r  =(x^+V^+z^)^  0  =  cos""^  — ^  =  sin~^  — ^ 
Tq  '■^0  ^0  0^  ’  0  r  r 

then  a  point  of  position  vector 

f.  H(p.,(().,Z.) 

with 

Pi=P0’  +  "i="0 


(2.2) 


(2.3) 


lies  on  Sj. 
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Points  related  like  and  will  be  called  "corresponding  points". 
Such  points  are  illustrated  in  figure  2.2.  Wherever  and  fj  or 

their  similarly  subscripted  coordinates  appear  in  the  same 
equation,  they  are  to  be  understood  as  corresponding  in  this  sense. 
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Figure  2.1:  The  ensemble  of  bodies  with  surfaces  Sq, 
disposed  with  angular  periodicity  about  the  z  axis.  A  plane 
is  incident  from  the  direction  of  spherical  coordinates  6\(t)‘. 
scattered  field  is  observed  from  the  direction 


1’ 

wave 

The 
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Figure  2.2:  "Corresponding"  points  and  r.  on  bodies  Sq  and  Sj. 

The  vectors  c^C^q)  and  are  rotated  relatively  to  one  another 

about  the  z  axis  in  the  same  way  as  r^  and  f.. 
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3.  The  incident  field  and  its  harmonics 

In  this  section  the  incident  electromagnetic  field  is  defined.  The 
manner  of  its  decomposition  into  a  Fourier  series  is  described  and 
the  rotation  dyadic  is  introduced. 

A  plane  wave  is  incident  on  the  ensemble  of  scatterers  from  the 
direction  9^(t)^  see  figure  2.1.  The  electric  field  vector  at  a  point 
(x,y,z)  or  (p,(t),z)  is 

•  •  jk^(xsin9^cos(t)^ +ysin0*sin(t)^ +  zcos9^) 

^ 

;  jkr,(psin9*cos((})-(t)^)  +  zcos0^) 

=  e  E^e  ^ 

(3.1) 

where  kQ  is  the  wave  number  for  e  ^  time  dependence,  the 

factor  ^  is  suppressed,  and  is  a  unit  vector  specifying  the 
polarisation.  In  this  work  we  consider  the  two  cases  e  =  9  and 
Any  other  may  be  resolved  into  components  in  these 

directions.  We  note 

9^  =  X  cos 0^  cos (j)^  +  y cos 0^  sin  (|)^  -  zsin  9^ 

=  p  cos  0^  cos((j)  -  (j)^ )  -  ^  cos  0^  sin((j)  -  (|)^ )  -  z  sin  0^ 

and 


=  -X  sin  (j)^  +  y  cos  (})^ 

=  p  sin((j)  -  (j)^ )  +  ^  cos((j)  -  (t)^ ) 


(3.3) 


in  which  the  unit  vectors  p  and  ^  are  evaluated  at  the  observation 
point  (p,(t),z).  Here  and  subsequently  the  superscript  i  specifies  the 
incident  field  while  the  subscript  i  is  an  index  specifying  the  body 

Si- 

(1),  (2)  and  (3)  give,  for  e^=0^ 
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— ;  :  ;  jk„psin9^  cosCd)  -  (1)^) 

E'  =  (p  cos  9^  cos((t)  -  (])*  )e  ^ 

i  ;  jk^psin9^  cos((l)-(b^) 

-(i)cos9‘sin(({)-(t)’)e  ^ 

;  jk„psin9^cos((l)-(l)^)  jkp.zcos9* 
-zsin9^EQe  ^  )EQe  ^ 


(3.4) 


and  for 


•  .X  xi^  jkopsin9^cos((l)-(l)^) 
E  =(psin((t)-(t)‘)e  ^ 


+  ^cos((|)-(l)^)e 


j  jk„psin9^cos((t)-(t)^)  jk„zcos9^ 

-rSMp  ^  EqC  ^ 


(3.5) 


The  products  of  trigonometric  and  complex  exponential  functions 
in  (3.4)  and  (3.5)  are  expanded,  using  the  formulas  in  Appendix  A, 

as  Fourier  series  in  (|)-(1)^  Thus 


E'(f)=  I  E;^(f) 
n  =  — oo 


trl  /z:\ 


(3.6) 


where 


n'-'  -9  n' 


E^„(f)  =  E.cVf)e 


(3.7) 


will  be  referred  to  as  the  "n  th  harmonic"  of  E^  and 
for  =9^, 


ci(r)=cl«(f) 


n  ^  i  ^  n  jk^zcos9^ 

=  j"(-pcos9^jJ  (P^)-(j)cos9^— r-J  (P^)-zsm9^J  (P^))e  ® 
n  d1  n  n 


(3.8) 


while  for 


i  id)  n  n  ;  -  ;  jk^zcosG^ 

5;i(?)  =  c|j^(f)  =  -f(p— ®  (3.9) 


in  which 
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P*  =  k;Qpsin0\ 


(3.10) 


and  the  non-constant  unit  vectors  p  and  ^  are  evaluated  at  r. 

Consider  corresponding  points  and  r.  and  the  vectors  cJ^(rQ)  and 
cj^(rp,  see  figure  2.2.  The  position  vectors  are  related  through 


r-fo-Ri 


(3.11) 


where  Ri  is  a  dyadic  which  rotates  a  vector  about  the  z  axis 
through  angle  Oj.  In  rectangular  coordinates  (3.11)  corresponds 

to  the  matrix  equation 


r.  cosO.  -sinO.  Ol 

IX  1  1  Ox 

r.  ,  =  sinO.  cosO.  0 

ly  1  1  Oy 

r.  0  0  1 

.  izj  ■-  -‘I  Oz 


(3.12) 


where  rj^^  etc  are  the  rectangular  components  of  r..  From  their 

definitions  in  (3.8)  and  (3.19)  it  is  clear  that  the  are  related 

n 

similarly: 


(3.13) 


and  with  (3.7)  that  the  harmonic  fields  are  related  by 


— i  — i  = 


(3.14) 


4.  Induced  surface  currents 


In  this  section  the  harmonic  surface  current  densities  are  defined 
and  the  relations  between  the  currents  on  the  different  bodies  are 
deduced. 

The  nth  harmonic  of  the  incident  field  excites  a  current  at  f.  on  Sj 

of  vector  surface  current  density  harmonic 

current.  The  total  current  is  the  sum  of  the  harmonics. 

Sj  is  rotated  relative  to  Sq  through  the  same  angle  O.  as  the  field 

ij^lrp  is  rotated  relative  to  Ej^(fQ),  see  (3.14).  Thus,  relative  to  its 
own  orientation,  Sj  experiences  the  same  incident  field  harmonic 

as  does  Sq,  except  for  the  phase  factor  e  f  Therefore  the  vector 
currents  exhibit  the  same  relations  as  the  fields: 


-c 

J-  (f.)  =  Jn  (fn)  Rie 
in^  1  On^  0 


(4.1) 


(4.1)  reveals  a  major  benefit  of  expressing  the  field  as  a  sum  of 
harmonics.  As  the  currents  on  the  surfaces  are  simply  related  it  is 
necessary  to  determine  the  current  on  only  one  of  them.  If  the 
current  is  found  by  the  moment  method  to  be  discussed  below, 
and  it  is  expressed  as  a  set  of  N  unknowns  per  surface,  a  direct 
approach  would  require  the  solution  of  NI  equations  in  NI 
unknowns,  the  solution  time  being  asymptotically  proportional  to 
(NI)^.  With  the  present  approach  the  solution  time  is  proportional 
to  per  harmonic.  We  shall  presently  see  that  the  number  of  - 
significant  harmonics  approximates  2k  times  the  ensemble 
diameter  measured  in  wavelengths. 
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5.  The  electric  field  integral  equation 

In  this  section  we  derive  the  electric  field  integral  equation,  the 
solution  of  which  is  the  surface  current  density  distributions  on 
the  bodies.  It  is  found  necessary  to  calculate  the  current  on  only 
one  body,  for  only  one  direction  of  incidence,  and  that  the  currents 
on  the  other  bodies,  for  other  directions  of  incidence,  may  be 
derived  from  the  first. 

The  currents  induced  on  the  surfaces  radiate  the  scattered  field. 
The  total  field,  incident  plus  scattered,  satisfies  the  boundary 
condition  that  everywhere  on  the  surfaces  its  tangential 
component  is  zero.  This  is  expressed  in  an  integral  equation 
whose  solution  is  the  surface  current  density  distribution. 

The  field  E^(f)  at  a  point  of  position  vector  r  radiated  by  the  nth 
harmonic  currents  on  all  the  bodies  is  the  sum' 

_  I-l  -  = 

E"„(r)=  I  JJJ?^(fj').Go(f,fj')ds'  ,  (5.1) 

i  =  0  f.'  onS. 

1  1 

where  ^(f,F)  is  the  dyadic  Green's  function  relating  the  field  at  f 
to  a  point  source  at  F . 

With  (3.11)  and  (4.1)  this  is  rewritten  in  the  form 

I“1  jnO.  =  ==  = 

E^(r)=  X  e'  ‘  (5.2) 

i  =  0  f^'cnSQ 

The  boundary  condition  is 

[E^(r.)  +  E^(r.)]  =0,  allr.  onS.,  j  =  0,l,...,I-l  (5.3) 

where  the  notation  [.]j  specifies  the  vector  component  tangential  to 


With  (3.7)  and  (5.2),  (5.3)  becomes 


I-l  jnO.  =  =  = 

[le  ‘  JJ  (Jl  (r„').Ri).Go(f|,rp'.Ri)ds']j 
i=0  P  J 


ro'onSo 


jn((j)j-<t)‘) 

=  [-E^cJ^Cf j)e  ^  ]^,  ail  Tj  on  Sj,  all  j 


(5.4) 


(5.4)  must  be  satisfied  for  all  j.  Explicitly  for  j  =  0, 


I-l  jnO.  =  =  = 

[  I  e  1  JJ(jS^(rQ').Ri).Go(fQ,rQ'.Ri)ds']^ 

i  =  0  fQ-onSQ 

•  jn((t)p, -(j)^) 


(5.5) 


With  (3.11)  and  (3.13),  this  may  be  written,  after  some 
rearrangement 


I-l  jnO.  .  _  _  _  _ 

[Ze  ^  J  IS  J^„(rg’).Ri  +  j,Rj‘.GoCfQ.fj').R.ds'], 

i  =  0  rg’onSg 


_  ;  W+j-'l'') 

=  [-E^c^(r.)e  ■'  ].,  r.onS. 

'■  0  ■‘t’  J  J 


(5.6) 


where  Rj  ^  is  the  dyadic  which  rotates  a  vector  through  the  angle 

-O;  and  in  which  it  has  been  noted  that  O. +0.  =  0.  ,  .  and 
J  i  J  i+j 


Using 


Rj  ^Go(rQ,rpRj  =  Go(r.,r.  _^  .' )  (5.7) 

(see  Appendix  B),  and  replacing  the  summation  index  i  by  i+j,  we 
obtain 


I-l  jn<I>.  =  = 

lie  1  JJ  J^^(rQ').Ri.Go(f.,rpds']^ 

i  =  0  f^'onS^ 


(5.8) 


(5.8)  is  the  boundary  condition  on  Sj,  Cf.  (5.4),  and  has  been 
derived  from  (5.5),  the  boundary  condition  on  the  body  Sq.  Thus 
satisfying  the  condition  on  Sq  satisfies  it  on  all. 


In  (5.5),  multiply  both  sides  by  and  change  the  order  of 

summation  and  integration  to  get 

"■q  (5.9) 

^  jn<!>n 


where 


=  I-l  jnO.  =  =  = 

G2^(ro,fQ')=  I  e  ^Ri.Go(rQ,rQ'.Ri) 


(5.10) 


and 


(5.11) 


The  azimuthal  direction  of  incidence  (|)^  does  not  appear  in  (5.9) 

and  hence  the  new  unknown  J^(fQ')  is  independent  of  (t)^  This  is  a 

further  benefit  of  the  present  approach:  if  we  are  interested,  as 
we  shall  be,  in  directions  of  incidence  over  a  range  of  azimuthal 

angles  (j)i,  all  with  the  same  polar  angle  0^  it  is  unnecessary  to 
solve  more  than  one  integral  equation,  (5.9).  The  currents  on  all 
bodies,  for  all  azimuthal  incidence  directions  are  then,  from  (4.1) 
and  (5.11), 
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Surface-fixed  coordinates 


Several  of  the  vector  quantities  of  concern  to  us,  such  as  the 

and  both  sides  of  the  integral  equation  to  be  solved,  (5.9),  are 
tangential  to  the  body  surfaces.  At  a  point  on  the  surface  the 
vector  may  be  specified  by  its  components  in  only  two  orthogonal 
directions.  In  this  section  these  properties  are  used  to  reduce  the 
apparently  three-dimensional  problem  to  a  two-dimensional  one. 

For  a  given  problem  we  define  a  pair  of  curvilinear  orthogonal 
coordinates  a, [3  over  the  surface  Sq.  This  serves  to  define,  at 

every  point  on  Sq  a  pair  of  orthogonal  unit  vectors  d(fQ),p{rQ) 
tangential  to  the  surface.  They  are  related  to  the  cylindrical  unit 
vectors  p,^,z  through  the  direction  cosines 

ttpCfo)  =  a(ro)-P(^0^ .  a^jj  =  a-i  >  “z  =  a-z .  Pp  =  P  p ,  Pjj,  =  P-i ,  P^  =  M 

(6.1) 

Note  that  all  the  quantities  in  (6.1),  excepting  z,  are  functions  of 
position,  as  shown  explicitly  in  the  first  relation. 

Each  side  of  (5.9)  is  decomposed  into  its  components  in  the  d  and 
P  directions; 


aao)-  « 

Jo'onSo 

PcIq)-  V 

rg’onSo 


(6.2) 


where 


nq  n  ' 
(6.3) 


■  iC  ic  iC 

c  ^  c  \  c  ^ 
np  nq  nz 


-iT 


^(i)  ^Z 


ri  =  a,p;^  =  e,(l), 


are  the  cylindrical  coordinate  components  of 

^  =  0,(}),  i.e.,  the  respective  coefficients  of  p,  z  in  (3.8)  or  (3.9),  Pp 

(T|  =  a,P)  etc  are  given  by  (6.1)  and  superscript  T  denotes  matrix 
transposition. 


1  5 


The  surface  current  density  J^(rQ')  is  also  expressed  in  a  and  p 
components  and  respectively  ; 


(6.4) 


The  a  and  p  components  are  related  to  the  cylindrical  components 
thus: 


(6.5) 


The  subscript  on  the  matrices,  here  and  elsewhere,  specifies  the 
arguments  of  the  functions  which  are  its  elements. 

The  integrand  in  (6.2)  may  now  be  written  in  component  form  as 
follows: 


"Gpp  Gp(ti  Gpzl  ttp  Pp 

=  G<i)p  G#  Gctiz 

[Ozp  Gz$  G^d(fp,fg')  [^3 

(6.6) 

The  elements  of  the  G  matrix,  Gpp  etc,  are  derived  in  Appendix  C. 

With  (6.1)  and  (6.3)  to  (6.6),  the  vector  equations  (6.2)  may  now 
be  written  in  terms  of  the  a,p  components: 


(6.7) 


where 


~°aa 

°P“ 


°ap 

°»-(r  f  ■) 


(6.8) 


and  as  before,  ^  =  0,(|). 


7.  The  moment  method  of  solution 
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The  current  density  distribution  is  found  from  the  solution  of  the 
integral  equation  (6.7).  In  this  section  we  describe  a  method  by 
which  this  equation  may  be  solved. 

Two  sets  of  basis  functions,  f^(rQ)  and  defined  over  Sq 

and  the  current  density  distributions  are  expanded  in  terms  of 
them: 


(ro')  = 


a  a 

1  , 
i  =  \ 


=  1 


a/^(ro 


') 


(7.1) 


For  an  exact  representation  of  the  currents  the  summations  must 
in  general  be  infinite.  Truncation  to  the  finite  lengths  Lp  thus 

implies  approximation. 


(6.7)  and  (7.1)  give 


I  a,g““(r„)+  i  b  g“P(rQ)  = 
£  =  1  '  '  ”  <  =  1  '  ' 


-E  (r  ) 


Mn  iC 

■^0® 


(7.2) 


^  =  0or(j) 


where 

g?’‘(ro)=  ’l  =  “.P;x  =  a.P  (7.3) 

Tq'  ohSq 

and  the  alternatives  for  ^  correspond  to  the  alternative 
polarisations  of  the  incident  field. 

The  integrand  of  (7.3)  is  composed  of  known  functions  and  hence 
the  integral  may  be  performed,  at  least  numerically.  With  all  the 

g^^(fo)  known,  (7.2)  is  now  a  pair  of  simultaneous  linear 

equations  in  the  +  Lp  unknowns  a^.b^.  The  number  of 

equations  can  be  made  equal  to  the  number  of  unknowns  by 
enforcing  (7.2)  in  (L^^  +  Lp)/2  independent  ways.  The  simplest 
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way  to  do  this  is  to  enforce  (7.2)  at  +  Lp)/2  different  values  of 
Tq,  a  technique  known  as  "point  matching".  A  more  general 
technique  is  to  multiply  both  sides  of  (7.2)  by  (Lq,  +  Lp)/2 
independent  functions  of  f^  successively,  in  each  case  integrating 
the  results  over  Sq.  The  use  of  these  weighted  integrals  gives  the 
technique  the  name  "method  of  moments".  Point  matching  is  an 
application  of  the  moment  method  in  which  the  weighting 
functions  are  impulses. 


A  possible  choice  of  basis  functions  is  the  set  of  rooftop  functions 
illustrated  in  figure  7.1.  The  surface  is  divided  by  the  a,p  grid  into 
approximately  rectangular  patches,  numbered  ^  =  1,2,...,L.  (We  are 
taking  =  Lp  =  L.)  The  ^th  function  for  the  current  component 

in  the  a  direction,  f^(fQ)»  varies  with  position  f^  in  the  patch  as  a 

triangle  with  unit  height  in  the  direction  of  flow  a,  and  as  a 
constant  in  the  transverse  (P)  direction.  The  function  is  zero 

outside  its  patch.  The  function  for  the  p  component,  f^(ro)’  is 

similar  but  varying  as  a  triangle  in  the  p  direction.  The  patches 
overlap  such  that  the  edge  of  each  is  the  centreline  of  a  neighbour 
except  at  an  edge  of  a  surface. 


The  overlaps  together  with  the  triangularity  ensure  that  the 
current  density  is  continuous  in  its  flow  direction.  Also,  for  a 
surface  of  zero  thickness  terminating  in  an  edge,  the  normal 
component  of  the  current  density  at  the  edge  is  zero.  In  the  case 
where  the  bodies  are  in  contact  and  must  be  separated  by  a  cut,  a 
special  basis  function  must  be  used  at  the  cut  to  ensure  current 
continuity  there. 


If  point  matching  is  used,  a  convenient  set  of  points  for  its 
application  is  the  set  of  points  r^j^;  k  =  l,2,...,L,  located  at  the 


centres  of  the  L  patches.  Enforcing  (7.2)  at  these  points  yields  the 
2L  simultaneous  linear  algebraic  equations 


L  Lr 

a  +  y  b 

L  Lr 

^  Pa  ^  K  PP 


:  -Ep.e  ’  (fp.,  ), 

0  na^  Ok'^ 


-E  (r  ) 

0  np^  Ok^’ 


(7.4) 


^  =  0  or  (j) 
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8.  If  the  bodies  are  in  contact 


In  this  case  there  is  only  one  body,  with  a  shape  exhibiting 
angular  periodicity.  In  this  section  is  described  the  special 
treatment  required  for  this  case. 


The  original  single  body  is  sectored  into  I  separate  bodies  by  cuts 
exhibiting  the  same  angular  periodicity  as  the  original,  see  figure 
8.1  The  currents  flow  on  the  surface  of  the  original  body;  there 
are  no  currents  on  the  surfaces  exposed  by  the  cutting. 

Across  the  lines  of  the  cuts  on  the  original  surface  the  normal 
component  of  the  current  density  is  continuous.  Provided  the  cut 
is  not  along  a  discontinuity  of  the  surface  geometry  such  as  a 
ridge,  the  parallel  component  will  be  continuous  also.  Thus  in 
general  we  take  the  total  vector  current  density  to  be  continuous 
across  the  cuts. 


Consider  the  cut  separating  bodies  0  and  1  and  the  current 
densities  at  points  close  to  the  edges  formed  by  it;  let  f-,  be  a 

0,62 

point  on  body  0  near  its  edge  2  and  L  be  on  body  1,  near  its 

edge  1,  see  figure  8.2.  The  two  edges  are  formed  by  the  cut.  In 
the  limit  as  the  points  approach  the  cut  and  one  another,  the 
continuity  remarked  in  the  previous  paragraph  is  expressed  by 


(8.1) 


Now  let  be  the  point  on  body  0  corresponding  to  f.  on  1. 
From  (4.1),  (5.11)  and  (8.1)  we  find 


Jn(^0,e2)  =  jn%e^)-V 


jnO. 


(8.2) 


-  in  words;  the  vector  current  densities  on  the  two  cut  edges  of 
body  0  differ  in  direction  by  angle  and  in  phase  by  n<I>i,  but 
are  otherwise  equal. 

This  property  of  the  currents  may  be  built  into  the  basis  functions 
so  that  the  solutions  exhibit  it  automatically.  For  example,  a  set  of 
rooftop  basis  functions  running  from  edge  1  to  edge  2,  see  figure 
8.3,  may  include  a  special  function  consisting  of  a  half-triangle  of 
strength  unity  at  edge  1,  and  a  second  half-triangle  of  strength 


jnO, 

e  ‘  at  edge  2.  The  current  described  by  an  expansion  in  these 
functions  will  have  the  required  scalar  behaviour.  The  required 

vector  rotation  will  be  present  if  the  surface  directions  d(rQ),  P(fQ) 
rotate  from  edge  to  edge  through  the  angle  d)  y  separating  the 
edges  in  azimuth. 

Thus  an  appropriate  choice  of  basis  functions  ensures  current 
continuity  across  the  cuts. 

It  is  remarked  that  the  progressive  phase  shifts  of  the  currents  in 
the  separate  bodies,  see  (5.11),  ensures  that  if  the  ensemble 
includes  the  z  axis,  the  currents  flowing  towards  the  axis  sum  to 
zero  at  the  axis,  as  required  by  Kirchhoffs  current  law. 


cut  surfaces 
(no  currents) 


(a)  (b) 

Figure  8.1.  (a)  "Ensemble"  of  contacting  bodies  separated  by  cuts, 

(b)  One  of  the  separated  bodies;  no  currents  flow  on  the  surfaces 
exposed  by  the  cuts. 


Figure  8.2.  Showing  current  densities  adjacent  to  cuts  and  their 
continuity  across  cuts. 


9.  The  scattered  field 


In  this  section  an  expression  is  found  for  the  field  scattered  from 
the  stationary  ensemble  of  bodies  for  given  directions  of  incidence 
and  scatter. 

The  scattered  field  at  a  distant  point  r^  =(r,0^,(t)^)  is  the  sum  of  the 
fields  radiated  by  all  the  current  harmonics  on  all  the  bodies;  thus 
from  standard  radiation  theory  [e.g.  7], 


eV)  =  Kq  X  I 


oo  I  —  1 

a  4(fp.(0'9'+$'$")e 

n  =  -ooi  =  Or.'onS. 

1  1 


(9.1) 


where 

jkQr 

and  Zq  is  the  intrinsic  impedance  of  free  space. 

With  (3.11)  and  (5.12),  (9.1)  may  be  written  in  the  form 


E^f^)  =  Ko  X  Xe 
n  =  — ooi  =  0 


~  I-l  jo(<t. -((>') 


r'onS; 


ds' 


(9.3) 


where 


A,  =  f^(VRi). 


n9 


=  9^(P(f  ’XR.),  A 


(9.4) 


By  use  of  (3.12)  and  the  relations  among  vector  components  in 
rectangular,  cylindrical  and  spherical  coordinate  systems,  see 
Appendix  D,  the  terms  in  (9.4)  may  be  manipulated  to  yield 


jk  A 


jkQpQ'sin0^cos\j/  jk^z^'ccsO^ 
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(9.5) 

and 


A  xC 

n(j) 


jVr.TS 


sin  ¥  + 4)0' cosy) 
jkQPQ'sin0^cos\|f  jk^z^'ccsB^ 


(9.6) 


where 


V  =  (<|)q’  +Oj-(})^) 


(9.7) 


and  Pq',^q',Zq'  are  unit  vectors  at  r^'  in  the  p,(^,z  directions 
respectively,  (zq'  is  of  course  z,  which  does  not  vary  with 
position.) 

(9.5)  and  (9.6)  contain  products  of  trigonometric  and  exponential 
functions  which  may  be  recognised  as  closely  similar  to  those  in 
(3.4)  and  (3.5).  Their  right  sides  are  similarly  expanded  as  Fourier 
series  to  give 


jk„A^ 

A  e  0  ^  ff  ' 


m  =  -00 


jni(4>o' +<!>;-$  ) 


^  =  0  or  (j) 


(9.8) 

where  are  respectively  similar  to  the  defined  in 

m  m  V  j  mm 

(3.8)  to  (3.9)  with  0\(t)^  replaced  by  0^,(1)^  and  f  replaced  by  Tq’ ; 

thus 


0 


JkQZ^'  COS0® 


(9.9) 
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•m/;  ,  rn 


3  jkoZo'cose^ 


(9.10) 


where 


=  k^On'sinG^ 


(9.11) 


It  will  be  convenient  below  to  express  the  vectors  in  terms 

of  their  a,j3  components.  These  are  (Cf.  (6.3)) 


“p 

Pp  PzJ  ll  ■ 

nz 


C=e,<|) 


(9.12) 


where  are  the  cylindrical  coordinate  components  of 

or 

i.e.,  the  respective  coefficients  of  p,  z  in  (9.9)  for  ^  =  0,  or 
(9.10)  for  ^  =  (j),  and  the  Op  etc  are  defined  in  (6.1). 


(9.3)  with  (9.8)  to  (9.12)  becomes  after  some  rearrangement 


E"(rS)  =  KQ  Z  I  e 
n  =  — oom  =  — oo 


— j(n(|)^ +  m({)^)  j(n  +  m)<I). 


fg'onSo 


(9.13) 


Finally  we  note  from  the  definition  of  <!).,  see  (2.1),  that  the 

summation  on  i  is  zero  unless  n  +  m  is  an  integer  multiple  of  I,  in 
which  case  it  equals  I.  Therefore,  set  m  =  kl  -  n  and  sum  over 
integer  values  of  k  to  get 
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E^(f^)  =  KQl  X  I  eJ 

k  =  — oon  =  — oo 


a  Z  (fn'  _  n  (^n’  )e'  +  cfT_  „  (?.'  )(|)^)e 


.:.se  .^SS,^S(1)  "^^o'ds’ 


r^'oriSo 


-nV‘o  ^•^''kl-n'‘0 


kl-n^'O 


(9.14) 
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10.  The  spectrum  of  the  scatter  from  a  rotating 
ensemble 

In  this  section  the  hitherto  stationary  ensemble  is  considered  to 
rotate  about  its  axis.  The  scattered  field  is  modulated  as  the 
scatterer  moves.  A  formula  is  found  for  the  spectrum  of  the 
modulated  field. 

Rotation  of  the  ensemble  of  bodies  about  the  z  axis  with  angular 
velocity  Q  in  the  ([)  direction  may  be  simulated  by  setting  and  (j)^ 
in  (9.14)  to  rotate  with  velocity  -Q  about  a  stationary  ensemble 
while  holding  their  difference  constant.  Thus  with  4)^  =  -Qt,  (9.14) 
becomes 


(10.1)  represents  a  spectrum  with  lines  disposed  on  each  side  of 
zero  at  intervals  of  IQ.  This  is  the  radian  frequency  with  which 
the  bodies  pass  a  fixed  point;  given  the  focus  of  this  study  on  the 
modulation  of  radar  scatter  by  the  blades  of  a  jet  engine 
compressor,  IQ  will  be  referred  to  as  the  "blade  rate".  Q  will  be 
called  the  "shaft  rate". 


IQ  is  also  the  frequency  with  which  the  ensemble  re-presents  the 
same  aspect  to  the  radar.  It  is  thus  the  fundamental  frequency  of 
the  modulation  of  the  scattered  field.  It  is  therefore  in  accord 
with  expectations  that  the  spectral  components  are  separated  by 
IQ. 
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(10.1)  represents  a  spectrum  centred  on  zero.  If  the  suppressed 

factor  e  ^  is  reintroduced,  the  spectrum  is  centred  on  the  carrier 
frequency  coq.  An  alternative  interpretation  of  (10.1)  is  that  it  is 
the  spectrum  of  the  received  signal  homodyned  to  base  band. 
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11.  The  spectrum  calculated  from  the  current  expansions 

For  an  exact  calculation  of  the  spectrum  the  currents  must  be 
known.  If  they  have  been  found  by  the  methods  described  in 
section  6  and  7  and  are  available  in  the  form  of  (7.1),  then  (10.2), 
with  (9.12),  may  be  written  as 


^kl,9  _K  I  y 

pS  “V  ^  ^ 


(11.1) 

All  the  functions  and  constants  in  (11.1)  are  calculable  using  the 
methods  described  in  previous  sections.  This  equation  is  therefore 
a  formal  solution  for  the  spectrum  of  the  electomagnetic  scatter. 
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12.  Properties  of  the  spectrum 

Sections  10  and  11  have  presented  formulas  from  which  the 
spectrum  of  the  scattered  Held  may  be  calculated.  However, 
certain  general  features  of  the  spectrum  may  be  deduced  without 
resort  to  exact  calculation.  These  will  be  described  in  the  present 
section. 

(10.1)  reveals  the  scattered  field  to  consist  of  a  spectrum  of 
frequency  components,  centred  on  zero  (or  cdq,  if  the  carrier 

frequency  is  included),  separated  by  intervals  of  IQ  (the  blade 
rate)  and  extending  indefinitely  in  the  positive  and  negative 
directions.  This  accords  with  expectations  for  the  periodically 
modulated  field. 

Although  the  spectrum  extends  indefinitely,  only  a  finite  band  is 
significant.  It  will  now  be  estimated. 


The  J^I^q')  the  solutions  of  (5.9).  They  will  be  small  when  the 

forcing  function,  the  right  side,  is  small.  From  (3.8)  to  (3.10),  with 
the  familiar  fact  that  the  Bessel  function  is  small  when  the  order 

exceeds  the  argument,  the  right  side  of  (5.9),  and  hence  j^(fQ'), 

will  be  small  when 


|nl>kQPQSin6^  (12.1) 

The  largest  value  of  Pq  is  p^ax  ’  maximum  radius  of  the 

2,k 

ensemble  of  bodies.  Using  this  in  (12.1)  and  replacing  Icq  by  — 

■0 

(where  Xq  is  the  wavelength  at  radian  frequency  coq)  we  have  for 
the  maximum  value  of  ln|,  say. 


l^^lmax 


2k 

— P 
X  '^max 


sin  9^ 


(12.2) 


Multiply  both  sides  by  Q  and  interpret  as  the  frequency 

separation  from  the  centre,  highest  and  lowest 

significant  current  harmonics.  Thus 


— p  QsinB^ 


(12.3) 


In  the  right  side  of  (12.3),  is  the  maximum  velocity  of 

any  point  on  the  bodies  towards  or  away  from  the  source  of 

illumination;  with  the  factor  —  it  is  recognised  as  the  maximum 

^0 

Doppler  shift  of  the  illumination  at  any  point  on  the  bodies. 

Thus  the  currents  induced  on  the  bodies  consist  of  a  spectrum  of 
frequency  components,  centred  on  the  illumination  frequency  coq, 
spaced  at  intervals  of  the  shaft  rate  Q.  and  extending  with 
significant  strength  to  the  maximum  positive  and  negative  Doppler 
shifts  associated  with  the  linear  velocities  of  the  bodies.  The 
spectrum  extends  indefinitely  beyond  the  Doppler  limits  but, 
through  the  behaviour  of  the  Bessel  functions,  the  strength 
decreases  rapidly. 

From  the  foregoing  it  is  concluded  that  the  integrand  of  (10.2),  due 
to  the  factor  j^(rQ')>  is  negligible  when  ]n|>|n|j^^  as  specified  in 
(12.2). 

Through  the  definitions  (9.9)  and  (9.10)  and  arguments  essentially 
the  same  as  the  foregoing,  it  is  found  that  the  dyadic  in  (9.14), 

negligible  when 

lkl-nl>^p  sine^  (12.3) 

'  '  1  '^max  ^  ' 

From  this  result  and  (12.1)  it  is  readily  found  that  the  maximum 
value  of  |k|  for  which  the  integrand  of  (9.14)  is  non-negligible  is 

Mmax=-5rPraax<™®‘+™®*) 

IAq 

Thus  k  =  ±|k|j^^  give  the  highest  and  lowest  terms  in  the 

summation  in  (10.1);  the  spectrum  of  the  scattered  field  extends 
above  and  below  the  centre  by 

|k|  Ii^  =  - — p  Q(sin0^ +sin9®) 

'  'max  T  '^max 

^0 


(12.5) 


which  is  recognised  as  the  two-way  Doppler  shift  associated  with 
the  highest  speed,  viewed  from  the  source  and  receiver,  of  any 
point  on  the  ensemble  of  bodies. 

|k|j^ax  number  of  lines  of  significant  strength  on  either  side 

of  the  centre  in  the  spectrum  of  the  scattered  field.  Under  some 
circumstances  (12.4)  delivers  a  less  than  one.  In  this  event 

the  scattered  field  is  little,  or  un-,  modulated.  Two  extreme  cases 
are  of  interest: 


(i)  0^=0^=O  (axial  backscatter): 

(12.4)  predicts  no  spectral  lines;  despite  the  approximations  used 
to  derive  (12.4),  this  result  is  rigorously  correct,  as  will  be  seen 
below,  provided  I  >  2.  Axial  backscatter,  when  the  ensemble 
comprises  more  than  two  bodies,  is  unmodulated. 


(ii)  0^=0^=-j  (incidence  and  scatter  in  the  plane  of  rotation): 
(12.4),  with  (2.1),  predicts 


|k| 


AttP  P  O, 

^max  _  ^max  1 


max  j 
=  d 


‘0 

max,^/2 


Xo/2 


(12.6) 


where  d  is  the  maximum  chord  distance  from  one  body  to 

max.  A,  /  2 

the  next,  measured  in  half-wavelengths. 


As  derived  here  (12.6)  applies  in  the  case  of  plane-wave  incidence 
and  scatter  in  the  plane  of  rotation.  Radar  illumination  of  jet 
engines  is  commonly  close  to  axially  incident  and  (12.6)  may 
appear  not  to  be  relevant.  However,  the  field  falling  on  the 
rotating  compressor  blades  after  diffraction  through  the  intake  is 
not  a  plane  wave.  It  may  be  decoihposed  into  a  spectrum  of  plane 
waves  which  includes  components  in  the  plane  of  rotation. 
Similarly  the  scattered  field  contains  components  in  the  plane  of 
rotation. 


Thus  (12.6)  is  a  significant  result  for  jet  engine  modulation 
generally:  the  number  of  spectral  lines  on  each  side  of  the  centre 
is  approximately  equal  to  the  tip-to-tip  blade  separation 
measured  in  half-wavelengths. 
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If  thsre  are  no  significant  spectral  lines 

and  the  scattered  field  is  unmodulated.  Here  the  bodies,  even  at 
maximum  radius,  are  separated  by  less  than  half  a  wavelength; 
the  ensemble  is  electrically  similar  to  a  continuous  disc. 


13.  Axial  backcatter 


This  is  the  case  where  e^=9^=0.  Without  loss  of  generality  the 
incident  polarisation  may  be  taken  to  be  0^ 

From  (3.8),  with  Appendix  A, 

f  =  _  i(0)  -  J„  +  ,(0»  - 1(«)  + 

(13.1) 


J  (0)  =  1  n  =  0 

n 

=  0  n  0 


and  hence 


(13.2) 


=0 

n 


n9i±l 


The  similarly  defined  quantities  see  (9.9)  and  (9.10), 

reduce  to 


^‘l=+j|(P±jW  ^ 

=0  m  ±1 

m  m 


(13.3) 


It  is  clear  from  (5.9)  and  (13.2)  that  is  non-zero  only  for 


n  =  ±l 


(13.4) 


which  is  therefore  a  condition  for  the  integrand  of  (9.14)  to  be 
non-zero.  With  (13.3),  a  further  condition  is  that 


kl-n  =  ±l 


(13.5) 


(13.4)  and  (13.5)  permit  the  integrand  of  (9.15)  to  be  non-zero 
only  with  the  following  combinations  of  I,n,k: 


3  6 


I  = 

1 

n 

=  -1 

k  = 

-2,  0 

n 

=  +1 

k  = 

0,  +2 

I  = 

2 

n 

=  -1 

k  = 

-1,  0 

n 

=  +1 

k  = 

0,  + 1 

I  = 

>2 

n 

=  +-1 

k  = 

0 

(13.6) 

From  this 

it  is 

seen 

that  for 

I  =  1 

and  I  =  2 

the  values  for  kl  are 

-2,  0,  +2.  The  spectrum  contains  components  at  (oq-2Q,  coq  and 

(Oq+2D. 

For  I  >  2  there  are  no  side  frequencies;  the  scattered  field  is 
monochromatic  at  coq.  The  axially  backscattered  field  is 
unmodulated  when  the  number  of  bodies  is  greater  than  two. 

These  properties  of  axial  backscatter  have  been  reported 
previously  [7]  in  a  different  context. 
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14.  Calculating  the  spectrum 

In  this  section  are  brought  together  the  steps  to  be  followed  in 
numerical  calculation  of  the  spectrum. 

1.  Define  a  set  of  rectangular  coordinates  with  z  axis  along  the 
axis  of  the  ensemble  of  bodies.  This  serves  to  define  the 
cylindrical  and  spherical  coordinates  also  in  use. 

2.  Choose  a  single  body  from  the  ensemble  and  call  its  surface  Sq. 
If  the  bodies  have  a  common  region  such  as  a  central  hub,  a  slice 
of  angular  width  2;t/I  must  be  cut  from  it  and  included  with  the 
chosen  body.  Sq  does  not  include  the  surface  exposed  by  the  cut. 

3.  Define  a  system  of  orthogonal  curvilinear  coordinates  a,p  on  Sq. 
This  serves  to  define  at  every  point  Fq  on  Sq  a  pair  of  orthogonal 

unit  vectors  dCr^XPlrQ),  tangential  to  the  surface. 

4.  At  every  point  Fq  on  Sq  determine  the  direction  cosines 

cCp  =  a-P  ;  «,{)  =  =  a-z ;  Pp  =  P-P :  =  P-l ;  P^.  =  P-z  ■ 


5.  At  every  point  on  Sq  determine  components 


^i0  ^i9  ^i0  ^i(t)  ^i<t) 

C  ,C,,C  .C‘.C,  C' 

np  n(j)  nz  np  n({)  nz 


see  (3.8)  to  (3.10). 

At  every  point  F^  on  Sq  determine  the  quantities 


C®  =-Ef.e^  +c^®a, +c^®a  ) 

na  0  ^  np  p  n<j)  (j)  nz  z’ 

C®  =-E  +c^®B  +c^®B  ) 


'np 


=-E^e'^  +c^^a, ) 

na  0  ^  np  p  n(t)  <j)  nz  z’ 

=-E  +c^^B  +c^^B  ) 

'^nP  0  ^  np^p  ^n(()^(|)  ^nz^z^ 


6. 


(These  are  needed  for  solutions  for  incident  polarisation  in  both 
the  6  and  (j)  directions.  If  only  the  first  (second)  is  needed,  only 
the  first  (second)  pair  of  quantities  need  be  calculated.) 

7.  For  every  pair  of  points  Tq^q'  on  Sq  compute  the  elements  of 


pp 

G  X 
P4) 

^pz 

<t)p 

(j)(l) 

G. 

(j)Z 

zp 

G  . 

z(p 

G 

zz 

(This  is  the  matrix  representation  of  the  dyadic  Green's  function 
Gj:,  see  (C19),  Appendix  C.) 

8.  For  every  compute  the  elements  of 

^aa  ^a|3 
Gq  Gqq 

L  pa  PP 

see  (6.8). 

9.  Choose  two  sets  of  basis  functions, 

f“(fQ),  f^(rQ);^  =  l,2,...,L 

to  represent  the  a  and  p  components,  respectively,  of  the  surface 
current  density  distributions.  (We  are  here  taking  equal  numbers 
of  the  basis  functions.  This  is  not  strictly  necessary.) 

10.  Choose  L  values  of  f^,  with  k  =  l,2,...,L.  Calculate  the 
elements  of  the  four  L  x  L  square  matrices 

in  which  hX  represent  the  combinations  aa,  ap,  pa,  pp  and  k  and  i 
take  values  1,2,...,L. 

11.  Assemble  the  quantities  calculated  under  headings  6  and  10  in 
the  matrix  equation 


and  solve  for  the  constants  a^,b^  ,  ^  =  . 

The  current  distribution  through  its  representation  by  a 
summation  of  basis  functions,  see  (7.1),  is  now  known. 

12.  Compute  the  integrals 

r^'onSo 

for  the  combinations  of  parameters 

Ti  =  a,|3;  %  =  0,(1);  |m (integer)] <--  -■. 

kQ/2 


13.  Compute  the  spectral  components: 


=  KqI  I  eJ 


n  =  -oo 


CX  /V  —  n  ry  j(kl  “ 

/"Vo  ff  ^  n,ot  ^0  ^o' 

(  I  JJ  f^(rQ  )  s<j,  e  ds 

^  =  1  fg'onSg  ['=kl-n,aj(j  .) 

T  "  S0 

P  B  ^kl  -  n  6  j(^  ~ 

+  X  b,_  Ji  fP(V)  30  ^ 

I  1  Tq  onSQ  l_^kl-n,pj^-^.^ 


for  all  |kl(integer)|< 
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If  there  is  interest  in  the  spectrum  beyond  the  Doppler  limits,  the 
computations  here  and  under  the  previous  heading  should 
continue  beyond  the  values  of  m  and  kl  indicated. 
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15.  Conclusions  and  proposals  for  further  work 
Conclusions: 

A  theory  has  been  developed  by  means  of  which  the  radar 
scatter  from  an  ensemble  of  bodies  with  angular  periodicity,  such 
as  a  compressor  stage  of  a  jet  engine,  may  be  calculated. 

The  theory  exploits  the  angular  periodicity  of  the  scatterer 
in  such  a  way  that  the  computational  load,  compared  with  an 
approach  which  does  not  do  this,  is  reduced  by  a  factor 
approximately  equal  to  the  number  of  bodies. 

The  theory  leads  to  a  solution  in  the  form  of  the  spectrum  of 
the  scattered  field. 

The  formal  solution,  without  resort  to  numerical  results, 
reveals  many  features  of  the  spectrum  such  as  line  separation, 
bandwidth,  number  of  lines  of  significant  strength  etc. 

A  numerical  approach  has  been  developed  in  sufficient 
detail  that  the  equations  might  be  solved  by  a  programmer  with 
little  or  no  familiarity  with  electromagnetics. 

Proposals  for  Further  Work: 

The  theory  described  here  contemplates  the  scatterer  as  an 
ensemble  of  bodies  in  free-space  illuminated  by  a  plane  wave. 

The  jet  engine  departs  from  this  ideal  model  in  a  number  of  ways: 
The  engine  compressor  stage  is  surrounded  by  an  engine  cowling. 
It  is  one  of  several  stages  close  enough  together  that  strong 
electromagnetic  interaction  is  to  be  expected.  The  stages  exist  in 
the  presence  of  stator  blades  which  themselves  constitute  periodic 
structures.  All  these  departures  from  the  ideal  are  expected  to 
have  their  effects  on  the  spectrum. 

It  is  proposed  to  continue  this  investigation  towards  a  better 
understanding  of  the  effects  mentioned  in  the  previous  paragraph. 
It  is  believed  that  the  present  theory  constitutes  a  basis  for  such  a 
continuation. 


n  =  -oo  P  " 

V(P>4^^n-l(P>  +  J„  +  l<P» 
J^„(P)  =  (-1)"J„(P) 


Appendix  B 


THIS  PAGE  LEFT  UNINTENTIONALLY  BLANK 
PLEASE  REFER  TO  AUTHOR 
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Appendix  C  The  Dyadic  Green's  Function  for  the  Periodic  Array, 
G^. 

By  definition  of  the  dyadic  Green's  function,  the  field  at  a  point  of 
position  vector  r  due  to  a  point  source  of  vector  current  moment  p 
at  position  F  is 

E(f)  =  p.G^(r,f).  (Cl) 

Express  E(r)  in  terms  of  its  components  in  the  orthogonal 
coordinate  system  {11^112,11^}  with  unit  vectors  evaluated 

at  f,  thus, 

and  p  in  terms  of  the  system  vectors 

evaluated  at  F,  thus. 


Now  let  p  have  unit  strength  and  orientations 
successively  to  get 

G^^,(f,F)  =  fi(r).E(f;C(F))  (C5) 

where  E(r;^'(F))  is  the  field  at  r  due  to  a  point  dipole  at  F  of  unit 

A 

moment  and  orientation  1^',  and  rj  =  11^,112,113; C  =  CpC2’C3- 


From  standard  radiation  theory  [e.g.  8]  the  field  of  the  current 
distribution  JCF)  niay  be  found  from  the  formula 


E(r)  =  111  (G'J(F)  +  RG"R.J(F))dv' 

(C6) 

where 

G'  (R)  =  - - ^(-1  -  jkR  +  k“R2)e"J^ 

jcoe  4;tR-^ 

(C7) 

G'(R)=  ^  ^.(3  +  j3kR  k^R^le'j^ 

47cR^ 

(C8) 

R  =  f-F 

(C9) 

and  the  integral  is  over  the  volume  in  which 

J(F)  is  non-zero. 

(C6)  may  be  written  in  dyadic  notation: 

E(f)  =  JJJ  J(F).G^(r,F)dv' 

(CIO) 

where 

^(r,F)  =  G'!  +  G''M 

(Cll) 

and  I  is  the  identity  dyadic. 

Let 

E  =  E'+E'' 

(C12) 

where  E'  and  E"  are  the  parts  associated  with  G'  and  G" 
respectively  in  (CIO),  (Cll).  Then  using  (C5)  with 
{111,112,113}  =  =  we  have,  for  the  pp  element  of  the 

matrix  representation  of  G'l  in  cylindrical  coordinates, 

0'pp  =  p.B(f;p')  =  p.p'G'  (C13) 

where  p  and  p’  are  unit  vectors  in  the  p  directions  at  r  and  F 
respectively.  Evaluating  the  dot  product  in  (C13)  gives 

G'  =  G'  cosCp  -  (])’ ) 

PP 


(C14) 
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The  other  elements  are  found  from  similar  considerations,  and  are 
set  out  as  follows: 


G'  I  s  G' 


COS((t)  -  (j)'  ) 

-sin((j)  -  0' ) 
0 


sin((j)  -({)')  0 

cos(({)  -  (j)' )  0 
0  1 


(C15) 


Through  similar  application  of  (C5)  the  matrix  representation  of 
G''R^  is  found  to  be 

G”RR  =  G'’x 

’(pcos((t)  -  (t)' )  -  p'  )(p  -  p’  cos((j)  -  (|)' ))  psin((t)  -  ([)'  )(p  -  p'  cos((t)  -  ))  (z  -  z'  )(p  -  p'  cos((})  -  (j)' ) 
(pcos((})-(t)')-p')p’sin((j)-(j)')  pp'sin^((j)-(j)')  (z  -  z' )p' sin((j)  -  (j)' ) 

(pcos((t)-({)')-p')(z-z')  p' sin((t)  -  (j)' )(z  -  z' )  (z-z')^ 

(C16) 

The  matrix  representation  of  G^Cr,?),  see  (Cll)  is  the  sum, 
element  by  element,  of  the  right  sides  of  (C15)  and  (C16). 

To  compute  G^  in  (5.11),  consider 

=  I-l  jnO.  _  _  _ 

p.Gj;(r,fQ’)=  le  1  p.Ri.Go(f,fQ'.Ri)  (C17) 

i  =  0 


which  is  the  field  at  r  due  to  the  periodic  circular  array  of  point 

_  jn<I>._  = 

dipoles  of  vector  moment  p-=e  ^p.Ri,  i=l,2,...I-l,  at  the 
corresponding  points  F.  =  FQ.Ri,  see  figure  Cl.  Note  that 


jnO. 

p.  =  e  1  (p(r.'  )pp  +  (t)(r.'  )p^  +  z(f )p^) 


(C18) 


in  which  it  has  been  recognised  that  the  unit  vectors 
p(fp,  ^(rp,  z(rp  are  related  to  their  values  at  through  the 

rotation  Ri  and  the  scalars  do  not  vary  around  the 

P  T  ^ 


array. 


Then  with  (C15)  and  (C16),  may  be  expressed  in  matrix  form 
thus 

r  ,  .  ^  r  C.  S.  O' 

_  I  - 1  jnO.  1  1 

Gj-Cf.rg)*  Xe-'  ‘{G'(Rj)-S;  Cj  0 

i  =  0  0  0  1 

(pC.-p’)(p-p'Cp  pSj(p-p'C.)  (z-z'Xp-p’C 
+  G”(Rj)  (pC.-p')p’S.  pp'S.2  (z-z')p’Sj 

(pC.-p')(z-z')  p'S.(z-z’)  (z-z-)^ 

(C19) 

where 

C.  =  cos((t)  -  S.  =  sm(({)  -  (j)p  (C20) 

and 

I 

R.  =  f-f.  =(p^+p.'^-2pp.'cos((t)-(t)-')  +  (z-z/)^)^  (C21) 

Since  P-'  =  Po’>  =  V"  V’  simplified  to 

R.  =(p^+pQ'^-2ppQ'cos((l)-(i)Q'-Op  +  (z-ZQ')^)2  (C22) 


Appendix  D  The  expressions  f.(A.R),  0.(A.R),  (|).(A.R) 

A  is  an  arbitrary  vector,  R  is  a  rotation  dyadic  of  the  type 
defined  in  (3.12),  and  (3.12)  with  rotation  angle  O,  and  f,  9,  ^  are 
unit  vectors  associated  with  the  spherical  coordinate  system, 
evaluated  at  the  point  (r,  0,  (j)).  The  unit  vectors  may  be  expressed 
in  terms  of  unit  vectors  in  the  rectangular  system  thus; 

f  =  xsin0cos(j>  +  ysin0sin(j)  +  zcos9 

0  =  xcos9cos(])  +  ycos0sin(})-zcos0  (Dl) 

A 

(j)  =  -xsin(j)  +  ycos(j) 

Performing  the  rotation  according  to  (3.12)  on  the  vector 

A  =  xA  +  yA  +  zA  (02) 

x  y  z 

and  then  forming  the  dot  products  with  f,  §,  $  as  given  in  (Dl) 
successively  gives 

f .  ( A.  R)  =  A^  sin  0  cos((l)  -  O)  +  A^  sin  0  sin({l)  -  0)  +  A^  cos  9 

0.  (A.  R)  =  A^  cos  9  cos((t)  -  O)  +  A^  cos  0  sin((t)  -  O)  -  A^  sin  0 

^.(A.R)  =  -  A  sin((l)  -  O)  +  A  cos(<])  -  O) 

X  y 

(D3) 

Applying  the  first  of  (D3)  with  f  =  =  (1,0^,(})^),  ^  = 

and  R  =  Ri  with  rotation  through  angle  O.,  gives,  after  minor 

manipulation  involving  trigonometric  identities  and  the  relations 
between  spherical  and  cylindrical  coordinates, 

(F  Q .  Ri )  =  p'  Q  sin  0^  cos((j)'  Q  +0 .  -  ({)^ )  +  z’  Q  cos  9^ )  (D4) 

Now  let  0  =  9^  and  A  =  I^(Fq),  but  write  the  rectangular 
components  of  A  in  terms  of  the  cylindrical  components  of  P, 
thus 


Making  these  substitutions  in  the  second  of  (D3)  gives 


0^.(J^.Ri)  =  cosB^ cos(({)'q  - (j)^)  -  cosG^ sin((t)'Q  +0.  - 

-  sin  9^ 
nz 


(D6) 


Finally,  setting  ^  =  and  A  =  J^(Fq)  in  the  third  of  (D3)  gives  the 
result 


.Ri)  =  sin((l)'Q  +0.  -({)")  +  cos((t)’Q  +0.  -  (t)^ ) 

(D7) 
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